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We present a study of a classical ferrimagnetic model on a square lattice in which the two inter- 
penetrating square sublattices have spins one-half and one. This model is relevant for understand- 
ing bimetallic molecular ferrimagnets that are currently being synthesized by several experimental 
groups. We perform exact ground-state calculations for the model and employ Monte Carlo and 
numerical transfer-matrix techniques to obtain the finite-temperature phase diagram for both the 
transition and compensation temperatures. When only nearest-neighbor interactions are included, 
our nonperturbative results indicate no compensation point or tricritical point at finite temperature, 
which contradicts earlier results obtained with mean-field analysis. 
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I. INTRODUCTION 



Stable, crystalline room-temperature magnets with 
spontaneous moments are currently the subject of a great 
, deal of interest because of their potential device applica- 
■ tions, such as thermomagnetic recording [Q . It is widely 
lO ' believed that ferrimagnetic ordering plays a fundamental 
On ■ role in these materials and the synthesis of new ferrimag- 
nets is at the moment an active field in material science. 

In a ferrimagnetic material two inequivalent moments 
interacting antiferromagnetically can achieve a sponta- 
neous magnetization at temperatures that are low com- 
^ ' pared with the strength of the interaction. At these low 
Q , temperatures, the inequivalent moments are antiparallel 
O ' but do not cancel [^j3|. This is particularly obvious in 
^ ^ the case of a linear chain, where the sum of the moments 
in each unit cell can result in a large moment for the 
chain. If adjacent chains can be positioned such that 
their moments are parallel, then a transition can occur 
at low temperatures to a state of 3-dimensional (3-D) 
ferrimagnetic order [Q. 

Important advances have been made by several groups 
in the synthesis of ferrimagnetic chains (|,^. However 
it is difficult to achieve high critical temperatures with 
one-dimensional materials. Consequently the discovery 
of bimetallic molecular materials with spontaneous mo- 
ments at temperatures as high as 43 K has directed 
the experimentalist toward the formation of 2-D and 3-D 
bimetallic lattices 

Synthesis of single-chain and double-chain ferrimag- 
nets is now becoming standard, and attempts to synthe- 
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size higher-dimensional polymeric ferrimagnets are start- 
ing to give very encouraging results. Some of the materi- 
als currently under investigation are 2-D organometallic 
ferrimagnets 2-D networks of the mixed- metal ma- 
terial {[P(Ph)4]_[MnCr(ox)3]}„ where Ph is phenyl and 
ox is oxalate 3-D ferrimagnets with critical tem- 
peratures up to 240 K and the recently developed 
amorphous T/(TCA^£')a;-y (solvent) with ordering temper- 
atures as high as 400 K 

The intense activity related with the synthesis of fer- 
rimagnetic materials requires a parallel effort in the the- 
oretical study of these materials. Mixed Ising systems 
provide good models to study ferrimagnetism. The mag- 
netic properties of these model systems have been ex- 
amined by high-temperatures series expansions and 
renormalization-group JTsf , mean-field Q], effective-field 
approaches ||l5|,|l6|. An exact solution of a mixed Ising 
system on a union-jack lattice (equivalent to one of the 
models studied in this paper) has recently been found for 
a low-dimensional manifold in the parameter space [ p7[ . 
In this work we study a classical model of a ferrimagnetic 
system: a mixed spin 1/2 and spin 1 system on a square 
lattice. We are interested particularly in the phase di- 
agram and in the location and characterization of the 
compensation point: the one temperature where the re- 
sultant magnetization vanishes below the critical point 
(type N in the Neel classification Q). The behavior at 
the compensation point is of technological significance 
since at this point only a small driving field is required 
to change the sign of the resultant magnetization. In a 
model for a ferrimagnetic thin film, it has been found 
that the coercivity diverges at the compensation point 
p8| . Preliminary results on our model have been pub- 
lished elsewhere |l9| ]. 

Experimental studies on recently synthesized com- 
pounds such as N(n-C„H2„+i)4Fe"Fe"^(C204)3 with 
n=3~5 have found critical temperatures between 35 and 
48 K, and some of these compounds exhibit a compensa- 
tion point near 30 K [^ . 

In Sec. H we present the Hamiltonian of the model 
and its ground states. We next briefly describe the non- 
perturbative techniques used to study the model: Monte 



Carlo (Sec. Ill) and numerical transfer-matrix calcula- 
tions (Sec. 0. In Sec. we discuss our results and 
finally we present the conclusions of our work in Sec. [v| . 



II. MODEL AND ITS GROUND STATES 

Our model consists of two interpenetrating square sub- 
lattices. One sublattice has spins a on the lattice sites, 
where a has two states, (t=±1. The spins a are spin 1/2, 
but we choose to put the factor of 1/2 into the interac- 
tion parameters. The sites of the other sublattice have 
spins S which can have three states, S'=±1,0. Each spin 
(J has only S spins as nearest neighbors and vice versa. 
The Hamiltonian includes all possible nearest and next- 
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nearest neighbor interactions and external fields. It is 
given by 

v. = — Ji ^ (TiSj — J2 ^ o-jS'l — J3 ^ SjSi 

<nn> <nn> <nnn> 

-Ji ^ Cr.j(Tfe - J5 ^ S'^Sf ^ ^ {SjSf + S'jSl) 

<nnn> <nnn> <nnn> 

i j j 

where the sums X)(nn) ^^'^ S{nnn) 
nearest-neighbor (nn) and next-nearest-neighbor (nnn) 
bonds, respectively. The sums J2i ^^'^ J2j ^'un over all 
sites of the a and S sublattices, respectively. Each J is 
an exchange interaction parameter, D is the crystal field, 
and Hi and H1/2 are the external fields, all in energy 
units. 

From now on, we will label a model by enumerating 
the parameters different from zero in the Hamiltonian. 
For example the J1-J2-D model is the model with all the 
parameters in the above Hamiltonian zero except Ji , J2 , 
and D. In all cases, we fix Ji to be <0, so the coupling 
between the nn spins is antiferromagnetic. 

In order to find the ground-state diagram for finite val- 
ues of the parameters, we use a 2x2 cell. With rotational 
symmetry taken into account, it has 2^3^/2=18 configu- 
rations. In Table | we show the 18 different configurations 
of the unit cell with their respective energies and degen- 
eracies. Which of these states is the actual ground state 
depends on the values of the parameters in the Hamilto- 
nian. Figs. 1^ through ^ show the ground-state diagram 
for different combinations of parameters. In each graph 
the ground-state configurations are indicated using the 
notation employed in Table ^ The boundaries between 
the regions are obtained by pairwise equating the ground- 
state energies. 



III. MONTE CARLO CALCULATIONS 

Standard importance sampling methods pl| were ap- 
plied to simulate the model described by the Hamilto- 
nian of Eq. ([^) on square lattices of LxL sites with pe- 
riodic boundary conditions. Most of the data were ob- 
tained with L=40, but we also present some results for 
L—IO, 16 and L=60. Configurations were generated by 
sequentially traversing the lattice and making single spin- 
flip attempts. The flips were accepted or rejected with 
standard heat-bath dynamics. We use the very long pe- 
riod, on the order of 2^^, random number generator KISS 



(for Keep It Simple, Stupid) |22|. Data were generated 
with 25000 Monte Carlo steps per site after discarding 
the first 2500 steps. The error bars were taken from the 
standard deviation of blocks of 500 measurements each. 
We define P—l/k-QT, and take the Boltzmann's constant 
fce^l. Our program calculates the internal energy 
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U = ^{-H), (2) 



the specific heat 



C ^ ^[{H') - {n)% (3) 



the sublattice magnetizations AIi and M2 defined as 

(4) 





M2 = T^{y^S,). (5) 



and the total magnetization AI—l/2{Mi+gM2), where 
the factor 1/2 gives the correct normalization for the 
whole lattice since Mi and M2 are normalized for the 
sublattice. Throughout this paper we take the (7-factor 
to be 5=1/2. We also measured the order parameters 



(6) 



and the susceptibilities associated with M, Mi, M2, and 
0± . The averages are taken over all generated configura- 
tions, the sums over i are over all the sites with a spins, 
and the sums over j are over all the sites with S spins. 
There are L^/2 terms in each sum. We verified that our 
results are in agreement with exact enumeration studies 
for L—2, and that the ground-state diagrams are repro- 
duced for different combinations of the parameters in the 
Hamiltonian. 

For an infinite lattice the order parameter 0+ would 
not be defined with the absolute value in Eq. ^, and would 
change sign at the compensation temperature Tcomp- 
However, for a finite lattice the absolute values are re- 
quired to keep the order parameters nonzero in the limit 
of a long measurement time. An efficient way to locate 
Tcomp using the Monte Carlo data is to find the cross- 
ing point between the absolute values of the sublattice 
magnetizations, i.e., 

|Mi(r,omp)| =<?|M2(Teo,„p)| (7) 

with the conditions 

sign(Mi(Tcomp)) = -sign(M2(Tcomp)) and Tcomp < 71: 

(8) 

These relations assure that 0+(Tcomp) as defined in 
Eq. is equal to zero. 
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IV. TRANSFER-MATRIX CALCULATIONS 



Traditional numerical transfer-matrix (TM) calcula- 
tions were performed as a second nonperturba- 
tive method to obtain finite-temperature phase dia- 
grams, critical exponents, and compensation tempera- 
tures. These results were compared with the Monte Carlo 
results, as well as with previous mean- field calculations. 

For a square lattice with different spins on each of the 
two square sublattices care should be taken to ensure 
that the TM is symmetric. A symmetric TM is pre- 
ferred, since it is much easier numerically to calculate 
its eigenvalues and eigenvectors. We used two different 
TM constructions, both of which give symmetric transfer 
matrices, as detailed below. The largest eigenvalues and 
eigenvectors of these symmetric transfer matrices were 
then calculated using the NAG subroutine F02FJE. This 
subroutine requires only multiplication of an arbitrary 
vector by the TM, and consequently it is not necessary 
to store the entire TM in memory. This allows us to use 
very large transfer matrices in our calculations. For both 
TM implementations, the lattice is wrapped on a torus 
of finite width and infinite extent, and periodic boundary 
conditions are imposed. However, the periodic boundary 
conditions are for a given column of spins, and lead to 
different boundary conditions in terms of the primitive 
lattice vectors since the torus of spins may not be along 
a direction given by a single primitive lattice vector. 

Two different implementations of the transfer matrix 
were utilized. These are called TMl and TM2, and the 
details of the construction of the transfer matrices is 
given in the appendix. 

The remainder of the equations in this section are 
stated for the TMl implementation with N, but would 
be equally valid for the TM2 implementation with TV re- 
placed by N. The inverse correlation length, is given 
by the ratio of the largest and next-largest eigenvalues of 
the TM as 



(9) 



The scaling form for ^ is [g^ 

Cat - NTitNy^), (10) 

where t — \(T — Tc)/Tc\ is the reduced temperature. At 
t=0 Eq. ( [lO[) allows one to calculate the finite-strip esti- 
mates for the critical temperature Tc as the temperature 
where the phenomcnological scaling relation 



N N- 



(11) 



holds. Differentiating the scaling relation Eq. (10), and 
evaluating it at the estimated value of Tc given by 
Eq. ( pT| ) gives the finite-strip estimate of the critical ex- 
ponent yT='^/v as 1^ 



+ ^ = ln[N/{N+l)] ■ 

The differentiation in Eq. (|l2|) was performed as a two- 
point finite difference. 

The TM calculation of the compensation point was 
only done in the implementation TM2. The standard 
method of calculating the magnetization by diagonaliz- 
ing a 2x2 matrix formed from the expectation values of 
the magnetization operator using the two-largest eigen- 
vectors was used ^J]. If M is the magnetization operator, 
the solution of the equation 



det 



(l|M|l)-m (1|M|2) 
(2|M|1} (2|M|2)-m 



= (13) 



give the magnetization. Here {i\ is the left (right) 
eigenvector associated with the i*'' largest eigenvalue of 
the transfer matrix. For a given N the compensation 
temperature is the temperature below the critical tem- 
perature where the magnetization is zero. 



V. RESULTS 



We first tested the mean- field predictions that the Ji- 
D model has a tricritical point and a range of D values 
with a compensation point. The ground state for this 
model corresponds to the line J2=0 in Fig. ^ Figure ^ 
shows the finite-temperature phase diagram for the Ji-D 
model as a function of D/\Ji\ (remember Ji——1) ob- 
tained with the Monte Carlo and transfer-matrix meth- 
ods, and we also show in the same graph the mean-field 
results ||l4|. The transfer-matrix (TMl) results were ob- 
tained using Nxoo and {N -f 1)xcxd lattices, with N=A, 
6, and 8. The Monte Carlo data for the critical tem- 
perature were obtained from the location of the specific 
heat maximum for lattices with L—10 and L—16. These 
results were confirmed with L=40, whose results are not 
shown in the figure. The finite-strip width estimates for 
yT calculated with TMl and TM2 are consistent with the 
Ising value 2/t=1 and are presented in Fig. ^. Indeed, as 
N 01 N increases the value for j/t approaches the Ising 
value for all Ji| < 4. These numerical transfer- matrix 
results strongly suggest that there is a multicritical point 
only at T=0 located at the point Ji|=4. There is no 
indication of a tricritical point at finite temperature. It is 
possible that a tricritical point is located at a much lower 
temperature that we could study, but even at our lowest 
value of Tc we do not see any indication of a tricritical 
point in the behavior of ut- An effective- field calcula- 
tion fill also has found that Tc=0 for L'/|Ji|>4. Our 
detailed TM study shows that Tc goes linearly to zero 
as rc/|Ji|=2.45(±0.11)(4-D/|Ji|) as D/\Ji\^4:. This 
linear relation holds extremely well between D/| Ji|w3.8 
and our lowest value of Z3/| Ji |=3.9995 where we get 
Tc/|Ji 1=0.001225. 
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A further study of the finite-temperature phase dia- 
gram for the J1 — J2—D model leads us to the conclusion 
that there is no compensation point for any range of pa- 
rameters of this model. Adding the field interactions {Hi 
and H1/2) does not seem to change this fact. Thus we 
conclude that a compensation point can not be induced 
in this model by nearest-neighbor interactions. 

Our next step was to include next-nearest-neighbor in- 
teractions. Adding the J3 interaction (between the S 
spins), we found that the finite-temperature phase dia- 
gram for the Ji — J^—D model also showed no evidence 
of a compensation temperature. 

Based on the above evidence, we concluded that a 
model with a compensation temperature different from 
zero must have, at least, interactions with non-zero pa- 
rameters Ji and J4 (between the a spins). A general 
study of the finite-temperature phase diagram of the 
Ji — Ji—D model shows that a compensation point ex- 
ists for a certain range of the J4 parameter in the region 
of Fig. I bounded by Ja/\Ji\ > and D/\Ji\ < 4. In 
Fig. B we show an example of the behavior of the mag- 
netization where the compensation and the critical point 
can be clearly observed. In Fig. |^ we show the critical 
and compensation temperatures plotted against J4 for a 
particular value of D. The compensation temperature 
does not exist until the J4 interaction takes some mini- 
mum value, after which it is almost independent of J4. 
This minimum value depends on D as is shown in Fig. ||, 
where we plot as a function of D the value of J4 at which 
a compensation point (Tcomp < ^c) appears. Fig. ^ indi- 
cates that only at D/Ji=A can there be a compensation 
point without the nnn interactions (J4=0). However, as 
one sees in Fig. ^, the compensation temperature seems 
to go to zero at this point, as does Tc- 

The critical temperature calculated using the numer- 
ical transfer-matrix method (TM2) for the J1 — J4—D 
model, is shown in Fig. |o|(a). The Monte Carlo estimates 
for Tc obtained from the maximum value of the specific 
heat are in excellent agreement with the ones obtained 
with the transfer matrix. To facilitate the reading of the 
critical temperature we present the curves for some val- 
ues of D in Fig. |lO|(b). For large values of J4 the critical 
temperature seems to be independent of D, but for small 
values of J4 the dependence on D is clearly observed in 
the graphs. In Fig. |l^(c) we show a detailed view of the 
behavior of the critical temperature for a particular value 
of D calculated with different transfer-matrix sizes. The 
finite size effects can be clearly appreciated. Also in this 
figure we show an exact result calculated from . It is 
important to emphasize that the exact solution is only 
known for the combination of parameters shown by the 
dashed line in the figure. Our numerical techniques cover 
the entire region of parameters. 

The compensation temperature for the Ji — J^—D 
model is shown in Fig. where the dotted line corre- 
sponds to Fig. |9[ As before, we choose to plot the numer- 
ical transfer-matrix results. The Monte Carlo results are 
in excellent agreement. Again the data strongly suggest 
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that only at D/Ji=A is there a compensation tempera- 
ture for J4=0, and that it seems to be zero. There is no 
compensation temperature at J^—Q for any other value 
of D. 

Small scale studies that included the other parameters 
of the Hamiltonian seem to indicate that for small values 
of the parameters J3, J5, and Jg there are only minor 
quantitative changes in the behavior of Tc and Tcomp- 
However, the external fields Hi and H1/2 play a more 
important role and the same seems to happen with the 
parameter J2. This is expected since for large values of 
Z?, J2 plays a role similar to an external field. 

VI. CONCLUSIONS 

We have applied two nonperturbative methods: Monte 
Carlo and numerical transfer-matrix calculations to 
study a mixed Ising system on a square-lattice. The 
model has two interpenetrating square sublattices, one 
with spins a—±l and the other with spins ^=±1,0. The 
Hamiltonian has all possible fields and nn and nnn in- 
teractions. In order to study the ferrimagnetic behavior 
of the model, we choose the coupling between nearest- 
neighbors to be antiferromagnetic. We calculated exactly 
the ground-state phase diagrams. Also, we have obtained 
the finite-temperature phase diagram and the critical and 
compensation temperatures for some interesting combi- 
nations of parameters. We found excellent agreement 
between the Monte Carlo and numerical transfer-matrix 
data. Our results show that a compensation point is in- 
duced by the presence of an interaction J4 between the 
spin-1/2 spins a (next-nearest neighbors in the lattice). 
A minimum strength of the nnn interaction J4 > for a 
compensation point to exist was found to depend on the 
other parameters of the Hamiltonian. We have demon- 
strated this in particular for the crystal field interaction 
D. We found that the system with only nn interac- 
tions does not have a compensation temperature except 
at the point where the crystal field takes its critical value, 
D/l Ji|=4, and the compensation temperature and criti- 
cal temperature seem to both be zero at this point. We 
also failed to find any evidence of a tricritical point at 
a finite temperature. Our nonperturbative results are in 
contradiction with mean-field studies for the Ji-D model, 
in which a tricritical point at finite temperatures and a 
range of D values with a compensation point were found 
1^,^. Thus we expect that there may be regions in 
some experimental two-dimensional ferrimagnets where 
compensation points may vanish when the couplings be- 
tween nn and nnn spins are changed, for example by the 
application of external pressures. 
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APPENDIX A: APPENDIX: CONSTRUCTION 
OF THE TRANSFER MATRICES 



This appendix presents the details of the construction 
of the two implementations, labeled TMl and TM2, of 
the transfer matrices that were used in the numerical 
calculations. We use the notation and methodology of 
Ref. nil to show the form of the TM and that the TM is 
symmetric. It is important to realize that in this notation 
curly brackets denote the matrix product introduced in 
Ref. H. 

The first implementation of the TM (TMl) consists 
of N spins (T=±l in the first column and N spins 
5=0, ±1 in the second column. This structure is iter- 
ated for an infinite number of columns. This only gives 
a symmetric matrix only if the nnn interactions are zero 
( J3=J4= J5= J6=0). The TM can be written as the sym- 
metric matrix dJ^2ADi A"'"Djy2 where the D matrices 
are diagonal matrices. The 2^x3^ matrix A, written 
for A^=4, is explicitly given by 



A = <^ 



Q 



Q 



Q 
Q 



Q 
Q 



(Al) 



where the 2x3 matrix Q, which takes into account the 
interactions with Ji and J2, is 



Q = 



/ exp[/3(JiH 
I exp [/3(-Ji 



exp [/3(-Ji 
exp [(3{Ji - 



- J2) 
^2)] 



(A2) 



In the normal fashion, the 2 x2 diagonal matrix T)i/2 
contains interactions of the field i?i/2, while the 3^x3^ 
diagonal matrix Di contains interactions of the fields D 
and Hi. 

In the second TM implementation (TM2) each column 
contains N spins a and N spins S. The spins are num- 
bered so that in one column a spin a is the first spin, and 
in the second column a spin S is the first spin. The TM 
has thej'orin^D^/'^BDB'^D^/^ where each of the matri- 
ces is 6^x6^. The matrix B, written for A^=3, has the 
form 



B 



(Q 



R 



s 








s 




R 








s 


Q 


s 








R 




R 








s 


Q 


s 








R 





(A3) 



where the 2x3 matrix Q is given by Eq. (|A2D . The 2x2 
matrix S takes into account nnn interactions between 
spins a (interaction J4), 



exp (pJi) exp {—PJ4] 
exp (— /3 J4) exp {PJ4) 



(A4) 
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and the 3x3 matrix R takes into account nnn interac- 
tions between spins S (interactions J3, J5, and Je), 



/exp[/3(J5 + J3 + Je) 
R = 1 

V exp [/3( J5 - J3)] 



1 exp[/3(J5- Ja)] \ 

1 exp [/?( J5 + J3 - Je)] / 
(A5) 



The 6^x6^ diagonal matrix contains interactions be- 
tween the spins and the fields as well as nn interactions 
of Ji and J2 within a column. For A'^=3 it is given by 



D 



Q 

H QT 



Q 

H QT 



Q 

H 



(A6) 



where the 2x2 diagonal matrix h has diagonal elements 
exp(/3i/i/2) and exp(— /3i7i/2) and the 3x3 diagonal ma- 
trix H has diagonal elements exp[/3(i?i — D)], 1, and 

-Hi — D)]. T he 6 ^x6^ diagonal matrix D has 



exp[, 



the same form as Eq. (A6), but starts with H rather than 
h as the first element in the matrix product. 
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FIG. 1. Ground-state diagram for the J1 — J2—D model. 
There are three regions, in each of which the configurations 
of the 2x2 cells are labeled as in Table |. The transition 
lines not parallel to the coordinate axes are labeled by the 
right-hand side of their defining equations, J2=aD—bJi 



FIG. 2. Ground-state diagram for the Ji — J3 — D model. 
The notation is analogous to Fig. 1. 

FIG. 3. Ground-state phase diagram for the J1-J4-D 
model. The notation is analogous to Fig. 1. 

FIG. 4. Finite-temperature phase diagram for the Ji—D 
model. The solid line corresponds to the mean-field approx- 
imation. The symbols O (*) with error estimates are 
Monte Carlo results on lattices with L=10 (16). The numer- 
ical transfer-matrix (TMl) results are also shown for system 
size iVxoo and (iV+l)xoo with iV=4 (Q), {+), and 

A''=8 (□). In the limit the model reduces to a nn 

spin ±1 Ising model, and the arrow indicates the exact value 
of Tc in this limit. 

FIG. 5. Numerical transfer-matrix results for the critical 
exponent yT~l/v. Finite-strip-width estimates using Nxoo 
and (iV-l-l)xoo lattices and TMl are shown for N=4 (Q), 
N=6 {+), and iV=8 (□). The resuhs obtained from TM2 
with Nxoo and (iV-l-l)xco are shown for N=2 (x) and N=3 
(★). It is clear that the results are consistent with the Ising 
value t/T=l. This becomes more evident as we increase the 
size. 
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FIG. 6. Magnetization vs temperature for the J1 — J4—D 
model at _D/|Ji|=3.6 and J4/|Ji|=3. The distinctive behav- 
ior of the magnetization at the compensation temperature 
and at the critical temperature can be observed clearly. The 
results shown are actually of | A/i | — l/2|Af2 1 for system sizes 
L=40 (+) and L=60 (x). The ones for M1 + I/2M2 and 0+ 
are almost identical except close to Tc where finite-size effects 
are largest. The solid curves give the magnetization from 
Eq. (hsf) from the numerical transfer matrix calculations with 
7V=2 through N=5. The numerical TM estimates for Tc from 
Eq. (^1]) are shown as vertical lines for N—2, 3 (dashed) and 
N= 3,4 (solid). In the insert we show the absolute value of 
the sublattice magnetizations, \Mi\ and 1/2 1 A/2 1. It is clear 
that at the compensation point the two sublattice magneti- 
zation cancel each other. In contrast, at Tc each one goes 
independently to zero, except for the remanent finite-size ef- 
fects. 

FIG. 7. Critical and compensation temperatures, Tc and 
Tcomp, for the Ji—Jn—D model at Z)/jJi|=3.6. The critical 
temperature is shown as a solid line and the compensation 
temperature as a dotted line. Monte Carlo results for L—40 
are shown by the symbol O with error bars. The numerical 
transfer-matrix (TM2) results for Tcomp with N=2 are repre- 
sented by X. For Tc the numerical TM sizes are N=2,3 (x) 
and A'^=3,4 (O)- The lines are guides for the eye. 

FIG. 8. Minimum value of J4 for the Ji — J4— D model for 
which a compensation point exists. The symbols have the 
same meaning as in Fig. ^. J4 (minimum) seems to go to zero 
as D/l Ji|^4. The lines are guides to the eye. 

FIG. 9. The compensation temperature at the values of 
J4 (minimum) given in Fig. |^, shown versus D. The numeri- 
cal transfer-matrix data, TM2 with N=2 (x), show that the 
compensation temperature approaches zero at the point given 
by J4=0, I>/|Ji|=4. 

FIG. 10. (a) Critical temperature for the Ji — Ji—D model. 
Tc is shown as a function of J4 and D. The values were 
obtained with the numerical transfer-matrix method TM2, 
with iV=2, 3. (b) Critical temperature for D/| Ji |=0, 1, 2, 3, 4 
and 5 taken from (a). The minimum of each graph depends 
on the values of J4 and D, but as IJ4I increases, Tc becomes 
independent of D. (c) A blow-up of the D/\Ji\=2 region 
is shown. This illustrates the finite-size effects going from 
A'^=2,3 (top solid line) to A''=3,4 (lower solid line). Also shown 
is the curve (dashed line) along which the exact solution is 
known, and the bicritical point (Q) fl^- 



FIG. 11. Compensation temperature for the J1 — J4—D 
model, calculated with TM2 and N=2. The dotted line cor- 
responds to Fig. y. 
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TABLE I. Ground-state configurations, degeneracies, and energies per lattice site for the 2x2 



cells. The spin states are indicated as follows: S=0 (0), S'=-|-l (©), S=—l (©), a=l (+) and 
a^^l (-). 



# 


Configuration 


Degeneracy 


Energy per site 


1 


+ 

+ 


1 


-El = —iHi — J4 


2 


+ 

— 


2 


E2 = J4 


3 


- 
- © 


1 


Ez = Hi/2/2 — J4 


4 


© + 
+ © 


2 


Ei = -i?i/2/2 -Hi/A + D/A -J1-J2-J4 


5 


© + 
— © 


4 


£5 = -Hi/ 4 + D/4 + J4 


6 


© - 
— © 


2 


Ee = Hi/2/2 - Hi/A + D/4 + Ji + J2 - J4 


7 


© + 
+ © 


2 


E7 = -^^1/2/2 + Hi/4 + D/4 + Ji - J2 - J4 


8 


© + 

— © 


4 


Es = Hi/4 + D/4 + J4, 


9 


© - 
— © 


2 


Eg = Hi/2/2 + Hi/4 + D/4 -J1 + J2-J4. 


10 


© + 
+ © 


2 


Eio = -Hi/2/2 + D/2 - 2J2 + J3-J4-J5 


11 


e + 


1 


Ell = D/2 - .h + Ji - .k 


12 


© - 
— e 


2 


E12 = Hi/2/2 + D/2 + 2J2 + J3-J4-J5 


13 


© + 
+ © 


1 


Ei3 = -Hi/2/2 - Hi/2 + D/2 - 2Ji - 2J2 -J3-J4-J5-J6 


14 


© + 
- © 


2 


Ei4 = -Hi/2 + D/2 -J3 + J4- - Je 


15 


© - 
- © 


1 


-E15 = -ffi/2/2 -H1/2 + D/2 + 2Ji + 2J2 -J3-J4-J5-J6 


16 


© + 
+ © 


1 


E16 = -Hi/2/2 + Hi/2 + D/2 + 2Ji - 2J2 -J3-J4-J5 + J6 


17 


© - 

+ © 


2 


En = Hi/2 + D/2 - J3 + J4 - J5 + Je 


18 




1 


Ei^ = H,_ ,/2 + Hi/2 ^ D/2 - 2.h + 2.1, - J-. - -h - + Je 
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